Ground state and excited states of a trapped dilute condensed Bose gas 
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I review recent theoretical treatments of a dilute interacting condensed Bose gas in a trap. Bogoliubov's classic 
results for a uniform condensate are generalized to include the effect of a trap, using the Gross-Pitaevskii 
formalism (for the condensate) and the Bogoliubov equations (for the linearized small-amplitude excitations 
of the condensate). Several recent theoretical studies are discussed along with some open questions.* 



1. INTRODUCTION 

Although 4 He was liquified in 1908 and the 
anomalous peak in the heat capacity at T\ s» 
2.2 K was discovered in 1932, its remarkable low- 
temperature superfluidity was observed only in 1938 
by Kapitza and by Allen and Misener jl], . Soon af- 
terward, Landau's bold 1941 theory of superfluidity 
of He II (including both the two- fluid hydrodynamics 
and the quasiparticle picture of phonons and rotons) 
|], ||, was brilliantly confirmed by Andronikashvili, 
Kapitza, and Peshkov jl], g]. In 1947, Bogoliubov 
Hi 01 proposed the first microscopic description 
of superfluidity in a uniform weakly interacting Bose 
gas and (with acknowledgment to Landau) suggested 
the correct generalization to a dilute gas with strong 
repulsive interactions. This work is summarized in 
Sec. 2; it has since served as the basis for studies of 
a nonuniform dilute Bose gas p] [| [To||. 

The last year has seen the exciting development 
of wholly new Bose condensates consisting of dilute 
alkali atoms @, @ @. At this conference, the 
F. London award has been presented to Cornell and 
Wieman, with comprehensive invited talks by Wie- 
man, Ketterle, and Jin. The present work seeks to 
build on these experimental achievements by review- 
ing the theoretical situation for a dilute Bose gas 
at low temperatures, well below the onset of Bose- 
Einstein condensation. The presence of the confining 
(typically harmonic) trap introduces essential new 
features that are summarized in Sees. 3 and 4 (for 
the ground state and low-lying excited states, respec- 
tively). Section 5 considers various specific calcula- 
tions that have already been performed, and Sec. 6 
provides a partial list of open interesting questions. 



2. REVIEW OF BOGOLIUBOV THEORY 

Five decades ago, Bogoliubov || [7] proposed a 
simple and intuitive description of a uniform dilute 
Bose gas, focusing particularly on how the interact- 
ing ground state differs from that of an ideal Bose gas 
jj|. Since the repulsive interactions favor a uniform 
configuration, Bogoliubov used periodic boundary 
conditions to ensure that the unperturbed ideal-gas 
single-particle ground state is the uniform Fourier 
component with k = 0. This choice is neces- 
sary because, for other boundary conditions, the un- 
perturbed ground state would not be even approx- 
imately uniform; in the thermodynamic limit of N 
particles in a volume V (keeping the particle den- 
sity n = N/V fixed), the introduction of weak re- 
pulsive interactions produces a dramatic change in 
any nonuniform equilibrium density profile (for ex- 
ample, an ideal Bose gas in a box with rigid walls 
has a density proportional to products of squared 
real trigonometric functions, whereas the interacting 
density is essentially uniform except very near the 
walls). 

An ideal Bose gas at zero temperature has all par- 
ticles in the condensate, so that the number No of 
particles with zero momentum is just the total N; for 
an interacting system, however, No is smaller than 
N (in lowest-order perturbation theory, two parti- 
cles can scatter out of the condensate and occupy 
the many zero-total-momentum states with separate 
momenta k and — k). For a dilute interacting Bose 
gas with a macroscopic condensate (No ^> 1), Bo- 
goliubov made the very simple observation that the 
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commutator of the second-quantized condensate op- 
erators [co,cJ] = 1 is small compared to their sepa- 
rate effect on the ground state (which yields a factor 
of order \/N~q ^> 1). Thus these condensate operators 
can be treated as simple numbers. In addition, Bo- 
goliubov assumed that most of the particles remain 
in the condensate, with N' = N — No <C N so that 
the total depletion of the condensate remains small. 
This condition means that it is sufficient to keep only 
terms of second order in the noncondensate oper- 
ators, and a linear (canonical) transformation can 
diagonalize the resulting approximate hamiltonian. 

For plane-wave states (oc e lk r ) at T = K, a 
particle with wave vector k has two characteristic 
energies: the kinetic energy 



(1) 



and the mean "Hartree" interaction energy with the 
macroscopic condensate 



V H = n V , 



(2) 



where Vq is the spatial integral of the interpar- 
ticle potential (namely the k = component of 
the Fourier transform); note that Vh is effectively 
a constant optical potential affecting the propaga- 
tion through the medium. As proposed by Bo- 
goliubov/Landau and proved subsequently by Lee, 
Huang, and Yang [|J [|, the correct generaliza- 
tion for a strong repulsive interparticle potential 
(whose Fourier transform diverges) is to replace Vo 
by Airah 2 /m, where a is the s-wave scattering length; 
in this case, the Hartree potential becomes Vh = 
Airah 2 no/m. 

One of Bogoliubov's central results is the energy 
eigenvalue for a plane wave: 



E k = x /2V H T k +T 2 . 



(3) 



Here, the first term under the square root is propor- 
tional to fc 2 and the second to fc 4 , so that Ek has the 
following limits 



E k 



^AnanQ h 2 k/m, 
h 2 k 2 /2r, 



2m, 



for k 
for k 



• 0: 



• oo. 



(4) 



In particular, the long-wavelength dispersion rela- 
tion is simply that of a sound wave (a phonon) 
with propagation speed s = ^/Anano Ti/m, given by 
the appropriate thermodynamic compressibility ||. 
This characteristic linear long-wavelength spectrum 
reflects the presence of the Bose condensate (note 



that s 2 oc no) and differs qualitatively from that of a 
dilute Fermi gas (for comparison, a nucleon moving 
through a nucleus with wave vector k has an energy 
« Vh + 7fe, which has a gap at k = 0). For a uni- 
form Bose gas, the transition between the phonon 
and free-particle spectrum occurs at a wavenumber 



fco = \/8irano — \ 2ms/Ti 



(5) 



obtained by setting T ko = Vh ■ The interactions must 
be repulsive (a > 0) to ensure that s and fco are real; 
in addition, the interactions play a crucial physical 
role, for fco vanishes if a — ► 0. As noted by Bogoli- 
ubov (based on Landau's quasiparticle picture), the 
critical velocity v c for the destruction of superfluidity 
is here simply the speed of sound, s oc •Jan®, again 
emphasizing the need for interparticle repulsion (in 
particular, v c = for a uniform ideal Bose gas). 

In addition to the energy eigenvalue, Bogoli- 
ubov's canonical transformation determines the cor- 
responding eigenstate, which involves a Bose "quasi- 
particle" created by the operator 



Vk C_k- 



(6) 



This description looks very similar to the familiar 
(and later) BCS theory of superconductivity |Q, but 
here the "coherence factors" satisfy the condition 



v 2 = l 



(7) 



to ensure that the quasiparticle operators obey Bose 
commutation relations. In the interacting ground 
state, Bogoliubov showed that the occupation num- 
ber N' k of the state with wave vector k is simply 
v\ = h[Eu 1 (VB + T k)-l], with the limiting forms 



K 



fc /\/8fc>l, for fc < fc , 
fco74fc 4 <l, forfc>fc . 



(8) 



The total fractional ground-state depletion of the 
condensate follows by summing over all nonzero 
plane waves: 



n 
n 



d 3 k , 



(9) 



It is important to recall that this finite zero- 
temperature depletion arises solely from the presence 
of the repulsive interactions. At finite temperature, 
thermal fluctuations induce additional depletion of 
the condensate |hJ, analogous to that for an ideal 
Bose gas. Equation (^|) shows that the Bogoliubov 
condition of small depletion (N' -C N) thus requires 



\/noa 3 <C 1; this depletion is nonanalytic in the scat- 
tering length a and cannot be obtained with any fi- 
nite order in perturbation theory (this is the prin- 
cipal advantage of the canonical transformation). A 
combination with Eq. (^) yields the alternative con- 
dition koa <C 1 for a dilute Bose condensate. 

An equivalent hydrodynamic description charac- 
terizes the same basic physics. In this case, the op- 
erator 

/4« V^(4 + c -k) (10) 

creates a density fluctuation with wave vector k; it 
has the same energy given in Eq. (||), even though 
it is a different linear combination than the quasi- 
particle operator in Eq. (|^). It is easy to show that 
these "phonons" are longitudinal and hence irrota- 
tional; as discussed below, vorticity occurs only for 
condensates containing discrete vortex lines. 

3. NONUNIFORM CONDENSATE 

Gross and Pitaevskii || extended Bogoli- 

ubov's theory to describe a rectilinear vortex line in 
an otherwise uniform dilute Bose gas. The conden- 
sate density no(r) is now spatially varying, and it 
convenient to introduce a "condensate wave func- 
tion" \&(r) whose squared absolute value is just the 
condensate density ^(r)! 2 = ^o(r). Evidently, the 
total number of condensed particles is given by No = 
J d 3 r |^(r)| 2 , thus fixing the normalization. In the 
present context of trapped alkali atoms, the Gross- 
Pitaevskii method must be generalized to include the 
trap potential V^(r); although the actual traps are 
anisotropic, the present treatment will consider only 
the simpler case of an isotropic harmonic trap with 



V(r) = |mwjr 2 



(11) 



In contrast to the situation for a uniform con- 
densate, there are now three separate contributions 
to the total energy: the kinetic energy 



T 



ft 2 V 2 
2m 



(12) 



the spatially varying Hartree energy of a single par- 
ticle with the nonuniform condensate 

Vff (r) = n (r)V = n (r), (13) 

m 

along with the trap potential V(r). The conden- 
sate wave function satisfies a nonlinear Schrodinger 
(Gross-Pitaevskii) equation 



{T + V + Vh)9 = h% 



(14) 



where fi is the chemical potential (given at T = K 
by fi = dE/dN). The nonlinearity arises from 
the Hartree potential, proportional to \^\ 2 ■ Equa- 
tion Jl4| ) is formally like the "Ginzburg-Pitaevskii" 
equation fill , but its interpretation is very differ- 
ent. Here, the Gross-Pitaevskii equation describes 
a nonuniform Bose condensate with small depletion 
at low temperature, whereas the Ginzburg-Pitaevskii 
equation is analogous to the Ginzburg-Landau the- 
ory of the second-order normal-superconducting 
phase transition jlTj and thus holds only near the 
transition temperature, where the (thermal) deple- 
tion is necessarily large. 

The trap potential introduces an additional 
length scale do associated with the zero-point mo- 
tion in the trap; setting htvo — h 2 jmd\ yields the 
familiar oscillator length 



do = ^Jh/muio. 



(15) 



For comparison, it is convenient to use the scattering 
length a and the "coherence" or "healing" length 



£ = 1/fco = 1/V87ran = h/V2i 



(16) 



that characterizes the distance over which the con- 
densate wave function heals back to its asymptotic 
value when subjected to a local perturbation (£ is 
analogous to the BCS or Ginzburg-Landau coherence 
length in a superconductor [|T7|). 

In the case of trapped 87 Rb fllj , these parameters 
are a ss 10 nm and do w a few /im, so that a -C d^; 



for the typical particle density n 



the 



coherence length is £ ss a few x 0.1 ^tm. If £ » do, 
then the healing length is large compared to the trap 
dimension and the system is nearly ideal; in this case, 
the selfconsistent condensate wave function is the 
Gaussian ground state of the harmonic oscillator 
(the corresponding condensate density no — \^g\ 2 is 
also Gaussian). 

For most experiments, however, the opposite con- 
dition a -C £ *C do applies, and the system is di- 
lute and interacting rather than ideal, so that the 
interparticle repulsion expands the condensate sig- 
nificantly; in this limit, the kinetic energy operator 
has only a small effect on Eq. (Q), leading to an 
algebraic equation 



{V + Vh-^tf = 0. 



(17) 



The resulting "Thomas-Fermi" (TF) approximation 
for the condensate density jl8], Q yields 



4nah 2 



rr> 



n TF(r) = V - V(r) 



(18) 



wherever the right-hand side is positive and zero else- 
where. In the present case of an isotropic harmonic 
potential, the TF condensate density has an inverted 
parabolic profile Uoxf(t) oc (R 2 — r 2 ) 6(R — r), where 
6 denotes the unit positive step function, R/do f» 
(lbNa/do) 1 / 5 1 characterizes the "condensate ra- 
dius" R, and \i ~ \1iu)qR 2 j d\ is the chemical poten- 
tial ^0|. It is often preferable to define the dimen- 
sionless ratio Na/do, which is large (small) for dilute 
interacting (ideal) Bose gas. As an aside, bulk liquid 
4 Hc is "dense," for it has £ s» a. 

4. SMALL-AMPLITUDE EXCITATIONS 
OF A NONUNIFORM CONDENSATE 

To treat the small-amplitude excited states in a 
trap, it is necessary to generalize the Bogoliubov 
quasiparticles to the case of a nonuniform conden- 
sate (this approach was originally developed for a 
vortex in a dilute Bose gas |J, but the inclusion 
of the trap potential is straightforward). Formally, 
introduce an operator (/) ( r ) that creates a noncon- 
densate particle at r [it is simply the fundamental 
second-quantized field operator V^( r ) with the con- 
densate operator and wave function ^(r) removed]. 
This particle operator can be expressed as a linear 
combination of quasiparticle operators 7J and j_j 
[compare Eq. (||)] , where the index j denotes one of 
a set of complete normal modes and —j denotes the 
time-reversed set of quantum numbers pi] . The cor- 
responding expansion coefficients are the Bogoliubov 
amplitudes (analogous to wave functions) uj(r) and 
Vj (r) that obey the coupled linear Bogoliubov eigen- 
value equations [0, ^2| 

(T + V + 2V H - fx)uj - V H Vj 



Vhu 1 



(T + V + 2V H - fi)vj = -EjVj, 



(19) 
(20) 



where the condensate density plays the role of a 
nonuniform potential through the Hartree contribu- 
tion Vjj- (The similar coupled equations for a super- 
conductor are usually known as the Bogoliubov-de 
Gennes equations [^3], Q.) The associated eigen- 
values Ej determine the quasiparticle energies in the 
presence of the nonuniform condensate, analogous to 
the Bogoliubov energies in Eq. (|^) for plane waves. 

An equivalent approach is to introduce hermitian 
small-amplitude hydrodynamic operators pq, Eg] 



$(r) 



p(r) = Vno(r) [0(r)+^t(r)], 
h 



[0(r)-0t( r) ] 



(21) 
(22) 



2miyno(r) 

for a density fluctuation and a velocity potential fluc- 
tuation at the point r [the associated velocity is sim- 



ply v(r) = V$(r)]. These hydrodynamic operators 
are merely linear combinations of the noncondensate 
field operators with known nonuniform coefficients, 
and the corresponding normal-mode amplitudes are 
similarly linear combinations of the Bogoliubov am- 
plitudes. In the present case of a stationary con- 
densate (the generalization to a condensate with a 
persistent current is straightforward pj|), the hydro- 
dynamic operators obey the coupled linear equations 



<9$ 
~dt 



Airah 



§ + V-(n V$)=0, 

^ Vn 

4m 2 no \f no 



(23) 



V 2 p 



(24) 

The first is the familiar linearized equation of parti- 
cle conservation, and the second is a linearization of 
the classical Bernoulli's theorem for a compressible 
irrotational isentropic fluid obtained from the appro- 
priate quantum energy-density functional [p5[ . The 
corresponding hydrodynamic amplitudes obey simi- 
lar coupled eigenvalue equations that determine the 
normal- mode frequencies ujj = Ej/h. 

The standard theory of linear response (27|] shows 
that the density-density correlation function charac- 
terizes the response of the condensate to a weak ex- 
ternal perturbation that couples to the density (ex- 
perimentally, this coupling can be obtained for a 
trapped Bose gas by modulating the curvature of the 
external trap p8| , p9|). Since the density operator is 
here merely a linear combination of the original par- 
ticle operators (this special form reflects the presence 
of a macroscopic Bose condensate, and the situation 
in a Fermi gas is very different), the resonant fre- 
quencies of the density-density correlation function 
are necessarily the same as those of the Bogoliubov 
equations, which provides an experimental method 
to determine the quasiparticle eigenvalues Ej. 

If treated exactly, the Bogoliubov and the hy- 
drodynamic descriptions both yield the same physi- 
cal information about the eigenvalues and eigenfunc- 
tions. The former relies on a quantum formalism 
with simple boundary conditions as r — ■+ 00, and the 
condensate density appears solely as a multiplicative 
coefficient. In contrast, the latter involves various 
spatial derivatives of the condensate density, but it 
offers the advantage of physical intuition (at least 
from a classical perspective), for the hydrodynamic 
amplitudes are like the normal modes of a selfgrav- 
itating nonuniform star. Indeed, the recent field of 
"helioseismology" uses the observed normal modes 
of the sun to infer the (otherwise invisible) behav- 



ior below the solar surface; similar studies may be 
applicable to the trapped alkali Bose condensates. 

5. RECENT THEORETICAL STUDIES 

The last year has seen the appearance of many 
theoretical papers and preprints on various aspects 
Bose condensation in dilute trapped gases |50|; the 
present section can only summarize a few selected 
studies that illustrate the recent trends. 

5.1. Stationary condensate 

The GP equation (|l4|) for a spherical static Bose 
condensate in an isotropic harmonic trap reduces 
to a nonlinear ordinary differential equation that 
can be solved in many ways, and the generaliza- 
tion to anisotropic traps is not difficult. A collab- 
orative group from Georgia Southern/NIST/Oxford 
[ 31 , 3^, [3^] has performed extensive numerical stud- 
ies of the radial condensate wave function (see also 
|34| ), including the dependence on the dimensionless 
parameter Na/do for both positive and negative val- 
ues (the latter corresponds to an attractive interac- 
tion). In addition, analytical studies have treated the 
same problem variationally f35| , [36f , using trial func- 
tions that interpolate smoothly between the (ideal) 
Gaussian limit (for Na/do *C 1) and the dilute (TF) 
limit (Na/d a > 1). 

Usually, the most relevant physical limit is 
Na/do 3> 1: when the TF approximation describes 
the spatial variation of the condensate density ex- 
cept in a thin layer near the surface, where the 
condensate kinetic energy becomes significant. Dal- 
fovo, Pitaevskii, and Stringari have recently studied 
this problem analytically using boundary-layer tech- 
niques |37j]. 

5.2. Rotating condensate (vortices) 

The existence and behavior of vortices in rotating 
superfluid 4 Hc has long been studied (see, for exam- 
ple, |38|, |3j|). In particular, the superfluid velocity 
v s (r) is proportional to the gradient of the phase of 
the condensate wave function. As a result, the su- 
perflow is irrotational "almost" everywhere except 
in discrete singular regions of dimension ~ £ (the 
nodes of \P) where the phase is undefined; superfluid 
vorticity (namely, the regions where V x v s 7^ 0) is 
localized in the cores of the superfluid vortices. In 
addition, W is single-valued whenever its coordinate 
r traverses a closed path in the fluid, which ensures 
the quantization of circulation k in units of h/m. 

The simplest system for liquid 4 He is a long cir- 
cular cylindrical container with radius R that can 
rotate about its symmetry axis. For small angular 



velocity f2, the fluid remains at rest with no vortices, 
but, at a critical value fi c i ~ (Ti/mR 2 ) ln(i£/£), it 
becomes favorable to create one singly quantized vor- 
tex line on the axis of the cylinder. Similar ques- 
tions have been considered for axisymmetric har- 
monic traps in the context of Bose condensed al- 
kali atoms p^, |33|, [54j, where, for current experi- 
ments on 87 Rb, the critical angular velocity should 
be f2 c i ~ 50 rad/s. 
5.3. Excited states 

As noted in the paragraph below Eq. ( |2J ) , the lin- 
earized normal modes of the condensate characterize 
both its small-amplitude dynamics and its linear re- 
sponse to a weak external perturbation that couples 
to the particle density. As a result, considerable the- 
oretical effort has been devoted to the calculation 
of these normal modes, both numerically |40[ and 
analytically in the TF limit of a large condensate 
pl[ . For an axisymmetric trap, the eigenstates can 
be labeled with the azimuthal quantum number m. 
In recent experiments [§8j, |2j|, the trap curvature 
was modulated at a prescribed external frequency, 
and the observed resonant behavior measures the ap- 
propriate eigenfrequencies. Detailed comparison be- 
tween theory and experiment yields good agreement 
for the lowest states with m = and m = 2. 

6. SOME OPEN QUESTIONS 

This field is developing rapidly, so that any se- 
lection of open questions is somewhat arbitrary (and 
sure soon to be irrelevant). Nevertheless, the follow- 
ing few topics are likely to remain of interest for at 
least the immediate future. 

6.1. Normal modes of a large condensate 

In the limit of a large condensate (Na/d 1), 
the TF approximation for the condensate density 
holds apart from the surface region. An incompress- 
ible spherical fluid has surface waves as its lowest 
normal modes, but it is less clear what happens for a 
dilute nonuniform Bose gas. In principle, numerical 
analysis can yield all the normal modes for a par- 
ticular configuration, but analytical results would be 
helpful to gain some insight into the relevance of var- 
ious physical parameters and the interplay between 
bulk modes and surface modes. 

6.2. Axisymmetric rotating condensate 

From an experimental viewpoint, the simplest 
way to create vortex lines in liquid He is to cool the 
sample while it rotates, so that the onset of the su- 
perfluid state occurs at fixed f2. This method works 
even for a circular cylindrical container, because the 



rotating walls entrain the viscous normal fluid above 
T\, thus creating the superfluid in a state of rotation. 

An alternative procedure is to cool a stationary 
container of He below T\ and then spin it up to 
some final angular velocity f2. For fl < fi c i, the su- 
perfluid typically remains at rest, but for > Q c i, 
the inevitable surface roughness of the rotating wall 
apparently initiates the formation of a vortex. Unfor- 
tunately, the detailed mechanism for vortex creation 
remains unclear, and the observed hysteresis implies 
significant metastability pH| . 

In liquid 4 He, only singly quantized vortices are 
expected theoretically, in agreement with many dif- 
ferent experimental observations p9[ ft2| . Imaging 
the vortices in rotating superfluid ^3e remained elu- 
sive for many years, in contrast to the corresponding 
success with quantized flux lines in type-II super- 
conductors [[17| . Nevertheless, clever use of ion trap- 
ping on the vortex cores (3^, |39| eventually yielded 
remarkably explicit images of the vortex positions 
43 [m]], which are really just the nodes in the con- 
densate wave function. Similar explicit images of 
vortices in dilute trapped alkali Bose condensates 
would be very desirable (it is not even obvious that 
only singly quantized vortex lines will occur). 
6.3. Nonaxisymmetric rotating condensate 

As noted above, spinning up a long circular cylin- 
der of superfluid He from rest frequently involves 
hysteresis, and it is not clear how a similar procedure 
for condensed alkali atoms in an axisymmetric trap 
can succeed in creating vortices (in the absence of 
physical boundaries, what does it mean to rotate the 
trap?). One obvious possibility is to use an asymmet- 
ric container (for example, different oscillator lengths 
d x and d y in the xy plane, with the harmonic trap 
rotating about the z axis). 

A similar question has been treated theoretically 
for an incompressible superfluid in a long cylinder 
with either elliptical or rectangular profile |45| . The 
principal new feature is that the rotating walls im- 
pel the superfluid into irrotational motion even in the 
absence of vortices. Indeed, for an elliptical profile 
with semiaxes a and b, the resulting induced angular 
momentum is reduced relative to that for solid-body 
rotation by a factor (a 2 - b 2 ) 2 /(a 2 + b 2 ) 2 p5| , which 
vanishes for a circular profile but is otherwise posi- 
tive and less than one. In addition, the critical an- 
gular velocity tt c i for the onset of vortex formation 
is increased relative to that for a circular cylinder for 
the following physical reason. It is easy to show that 
solid-body rotation $7 x r is the true equilibrium of 
an unconstrained rotating incompressible fluid in a 



container, but such a configuration is forbidden to 
a superfluid because v s must remain irrotational al- 
most everywhere (recall that solid-body rotation has 
uniform nonzero vorticity 2S~2). A superfluid in a 
rotating circular cylinder has no purely irrotational 
motion and can mimic solid-body rotation only by 
creating singly quantized vortices; in contrast, a ro- 
tating elliptic cylinder does have irrotational flow 
that already acts somewhat like solid-body rotation, 
reducing the need to create one or more discrete vor- 
tices. Careful experiments confirm this prediction in 
considerable detail |4(| . 
6.4. Excited states of a vortex 

A classical rectilinear vortex with circulation k 
in an incompressible fluid has oscillatory normal 
modes that propagate along (and are confined to) 
the axis, with a long- wavelength dispersion rela- 
tion us (Kfc 2 /47r) ln(l//c^), where k is the axial 
wavenumber and £ is the core size (assuming k^ <C 1) 
pSj , p9| . Similar vortex waves occur for a vortex in 
an unbounded dilute Bose condensate [l?J, where 
k = h/m. In addition, the phonon-like modes of a 
uniform Bose gas are modified by the presence of a 
vortex, leading to asymptotic phase shifts. The re- 
sulting phonon-vortex scattering contributes to the 
mutual friction, as observed in rotating superfluid 
4 He through an excess attenuation of second sound 
relative to that for stationary superfluid |3^, |3^| . 

The situation for one or more vortices in a 
trapped dilute Bose condensate will necessarily dif- 
fer considerably from that in an unbounded medium. 
As one simple example, the ratio of the condensate 
radius R to the vortex core £ is typically ~ 10, so 
that only relatively few vortices can form before the 
cores overlap and the system becomes "normal" (a 
finite-size analog of the destruction of superconduc- 
tivity at the upper critical field H C 2 [0). A funda- 
mental question concerns the various possible normal 
modes that can occur for vortex lines in these rela- 
tively small dilute trapped systems. Although there 
will certainly be phonon-like modes that are scat- 
tered by the vortices, it is less obvious that such a 
system can support an analog of the bound vortex- 
wave modes. 
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